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Abstract 

In this Letter we investigate the basic assumptions of the decay- 
ing dark matter (DDM) theory in the hght of recent advances in ob- 
servational and theoretical cosmology and physics, i.e. detection of 
massive astrophysical compact halo objects (MACHOs) and Super- 
Kamiokande results. Specifically, the consequences pertaining to the 
shape of the Milky Way galaxy dark halo are discussed. We find that, 
by taking into account the values of the main constituent of the mass 
in DDM theory, massive neutrino, with the mass of 30 eV, and life- 
time 10^^ s, the initially proposed value of extreme halo flattening 
q ~ 0.2 is no longer necessary, and that one can easily accommodate 
a much larger value of g ~ 0.6, that is in accord with all available 
observational data. 
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1 Introduction 



It is well known that every spiral galaxy consists of three main parts: central 
bulge, disk and dark halo (e.g. 0]). In this Letter we wish to investigate in 
more detail the idea that the halo of our Galaxy, Milky Way (in the further 
text also denoted as the Galaxy) is made mainly of two components: baryonic 
(MACHOs, see the discussion later) and non-haryonic (massive neutrinos 
with the mass my ~ 30 eV), and impact of such a composition of matter on 
the shape of the halo. 

The problem of the mass composition in spiral galaxies is well known and 
still unsolved (e.g. Ref. |]^). Rotation curve of the Galaxy (Ref. [Q) and 
non-existence of the Keplerian fall-off points out that there exists a large 
amount of matter that still has to be accounted for. It is a general opinion 
that a bulk of this matter resides in the dark halo of each spiral galaxy 
(Ref. [0, ^ |I[), although there are opposite standpoints that suggest the 
modification of the Newtonian dynamics (MOND) (Refs. ||21, ^ |23|) that 
recently produced quite good agreements with the observations of 15 rotation 
curves (Ref. [^). 

Massive neutrinos remain viable candidates for the "missing light" prob- 



lem (as defined by Schramm and Steigman ||3T|). It can be shown that under 
certain assumptions (Ref. ||41[]) massive neutrinos with the mass of the order 
of 30 eV can account for the mass in spiral galaxies (e.g. Ref. ||2^). These 
neutrinos play the central role in the decaying dark matter (DDM) theory 



usually associated with the name of Dennis W. Sciama |34 



One of the major motivations for undertaking this research is recent work 
of the present authors (Ref. |^), who marshalled lots of evidence pointing 
to a moderately flattened dark halo q ~ 0.6 for both our Galaxy and spi- 
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ral galaxies in general. Most of those arguments apply to an essentially 
baryonic (MACHO or gaseous) halo, but some are applicable to the general 
dynamically dominant halo component. Sciama's DDM theory is among the 
plausible and interesting theories predicting an extreme 0.2 flattening for 
the dark neutrino component. 

It should be noted that the presence of MACHOs detected through mi- 
crolensing searches does not automatically discard the DDM theory, neither 
its applications in the galactic context. It is still possible (albeit more and 
more difficult to maintain) that MACHOs are dynamically insignificant or 
only marginally significant. The estimates of the MACHO contribution to 
the cosmological mass density are still dependent on several not completely 
watertight assumptions [0, We shall discuss some possibilities for a 
combined DDM + MACHO model in the Section 4. 

2 DDM theory: Fundamentals 

According to the DDM theory it is assumed that neutrinos do have masses 
and that the more massive neutrino, ui radiatively decays into a photon and 
less massive neutrino, 1^2: 

Z/l ^ 7 + Z/2 (1) 

(e.g. Ref. B). If one further wishes to discuss the kinematics of this equation, 



one would obtain the following simple relation [Q: 

1 / ml \ 

where Greek letters in indices are related to the various types of neutrinos, 
i.e. tau, electron and mu neutrinos. It is also assumed that 

"^.e,. < 5 eV (3) 



therefore giving the following simple relation: 



^7 ~ (4) 



The basic values of the DDM theory are the following (Ref . |^ , updates 
in Ref. 



• According to the latest estimates the mass is: niy^ = 27 A ± 0.2 eV. 

• These neutrinos decay radiatively with the lifetime r = 2 ± 1 x 10^^ s. 
For the latest correction of this value, see the discussion later (based 
upon Ref. pi). 



• A decay photon energy is 13.7 ± 0.1 eV (obtained from the eq. [^]), so 
these photons can ionize hydrogen, but not helium. 

3 Flattening in the DDM theory and obser- 
vational data 

The significant fiattening of the halo was introduced by Binney, May and 
Ostriker (Ref. P|) who found that it is likely that if the kinematics of the 
objects in the halo are similar to the kinematics of the extreme Population 
II objects, then the massive halo should have the axis ratio c/a^O.5. One 
can now define the fiattening parameter ip as cosip = q = c/a, i.e. its cosine 
determines the shape of the halo En. The En notation is related to q as 
g = 1 - n/10. 

As for the DDM theory, Sciama proposed that the dark halo had to be 



extremely fiattened q ~ 0.2 — 0.3, in order to obtain scale height of electrons 
responsible for pulsar dispersion that is determined by the scale height of the 
dark matter (i.e. massive neutrinos). By scale height we assume the column 
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density on one side of the galactic plane divided by the volume density in the 
plane. Later on, Sciama [^] argued that the scale height of the ionized gas 
should be reduced to 900 ± 100 pc, thus eliminating the need of the strongly 
flattened halo. However, Nordgren, Cordes and Terzian found that the 
electron scale height is GTO^^Jq pc, and the mean electron density for the 
interstellar medium He = 0.033 ±0.002 cm~^. This is also in accordance with 
earlier results of Cordes et al. ^j. 

This leads us to important questions concerning the ionization problem 
in the spiral galaxies (e.g. Refs. 0, ^). We just mention here that the 
conventional sources such as O stars or supernovae have much smaller scale 
heights (~ 100 pc) than the aforementioned ones of the free electron compo- 



nent of the interstellar matter. Following Ref. we start with the equation 
that represents ionization equilibrium: 

— — = arig, (5) 

where n,y(0) is the neutrino density near the Sun, is the free electron 
density, r is the predicted value of the decaying neutrino lifetime and a is 
the hydrogen recombination coefficient excluding transition directly to the 



ground state (cf. Ref. [^). If one adopts a = 2.6 x 10 cm^s ^ and 
He = 0.033 cm~^ this leads to the following neutrino density: 

n^(0) = 2.83 X IOV23 cm"^ (6) 

where r is expressed as r = 10^'^r23 s. This can give the following value of 
the mass density Pu{0), assuming the neutrino mass of 27.4 ± 0.2 eV: 

p,(0) = 1.384 X 10-'V23 g cm-3 = O.O2r23A^0 pc'^ (7) 

where we, as usual, expressed mass density in the Solar masses (A^© = 
1.989 X 10^^ g) over cubic parsec (1 pc = 3.0857 x 10^^ cm). 
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Now we can see why the extremely flattened halo was needed in the DDM 
theory; namely using equation (0) one has: 



, n„\ 2 



(8) 



that gives, after inserting appropriate values, the value for the electron den- 
sity: 

0.016 < rie < 0.028 cm'^ (9) 

where we put Sciama's value = 2 x 10^ cm~^ (cf. Ref. [|3^) and let 
lifetime to take the values r = 2 ± 1 x 10^^ s. One can see that there are 
two ways to achieve agreement of theory with observations: one can either 
increase rii, as suggested by Sciama (Ref. by reducing the assumed scale 
height of the neutrino distribution of 8 kpc by a factor of ~ 4 to ~ 2 kpc 
which means that one must introduce extremely flattened halo, or, as we will 
show one can assume modified form of the equation (^) (see the discussion 
later). We first note that we obtained lower value than Sciama (Ref. ||36[| ) 
for Pv(0), O.O2r23A^0 pc"'^ vs. 0.03x23^^0 pc~^, because we used Nordgren 
et al. (Ref. ^6[) lower value for (0.033 cm~'^), rather than Sciama (0.04 
cm~^). In the Section, we shall see that this has profound consequences for 
flattening, even if we accept all other Sciama's premises. 

4 Uncertainties inherent in the flattening pa- 
rameter determination 

We would like now to investigate the general applicability of the decay- 
ionization equilibrium equation, and in order to do it, we would like to iso- 
late all possible sources of uncertainty. Therefore we propose the following 
improvements concerning equation (O): 



1) Additional ionizing sources should be taken into account, regardless 
on the assumed scale height. These can be O and B stars (the question of 
galactic disks opacity not being completely clear at present), as discussed, 
for instance, by Mathis |]l8l and others metagalactic background, cosmic 



rays, large-scale shocks, etc. 

In this case, the equation for decay-ionization equilibrium @ should read 

'^ + F = anl (10) 
r 

where F is the ionizing contribution of all other sources. Thus, we see that 

n,{<d) = {anl- F)t, (11) 

i.e. required neutrino density is decreased. Accordingly, pjy(O) is also de- 
creased, and necessity for flattening is decreased, as we shall see in more 
detail below. 

2) Assumption of "maximal neutrino" halo is unwarranted, even if de- 
sirable. The discovery of MACHOs, and substantial mass contribution of 
these objects to the required dynamical mass of the Galaxy (e.g. Ref. ||12|), 
immediately invalidates this assumption. So far, using the notation of the 
Sciama we have to use the decomposition 



Erot = Srot('^) + Erot(other) (12) 
Accordingly, the required lengthscale is decreased by a factor 

/ < 1. (13) 



Srot(z^) 



^rot 

This offers an interesting opportunity of building of more complex model 
incorporating other observed phenomena. If we suppose that the rest of 
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dynamical mass of the Milky Way halo is made of baryonic MACHOs of the 
same population as the one detected in the microlensing searches, then one 
should be able, in principle, to tightly constraint the parameter /, i.e. the 
decaying neutrino fraction of the dynamical mass. This is to be achieved 
by simultaneous consideration of (i) microlensing optical depths and their 
ratios in various directions, which would put constraints on the shape of the 
MACHO component ; (ii) ionization of Reynolds' layer as in the present 
paper; and (iii) constraints on the MACHO abundance stemming from the 
Big Bang nucleosynthesis constraints. 

We consider this to be a good bargain, since for price of introducing one 
additional parameter (/ ratio) in the theory, we get multiple advantages: 
accounting for observed microlensing events and baryonic dark matter on 
one hand, and non-baryonic dark matter on the other, as well as retaining 
some advantages of the classical DDM theory, like explanation of cosmological 
reionization. Therefore, in the course of future work, we shall try to show 
how this extension can be specifically realized. 

3) Column density uncertainty may vary within a factor or ~ 2) is 
introduced. 

4) Galactocentric distance of the Sun, as well as the rotational velocity 
of the LSR are still subject to some uncertainty (e.g. Ref.|Q). 

5) The least significant source of uncertainty is the uncertainty in the 
mass of decaying neutrino. It is manifested in transition from the neutrino 
number-density (as a discrete value quantifying ionization equilibrium) to 
the mass-density distribution. 

Now we pass to determination of the flattening parameter of the neutrino 



halo, using the method and notation of Sciama p6[. The vertical lengthscale 
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Figure 1: The flattening parameter g as a function of the electron density of 
the Reynolds' layer in the DDM theory. Vertical line denotes the recently 
obtained observational upper limit on the electron density. 
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can be written as 



[ani — F)Tmy 
Flattening parameter q can be approximated as 



Zq 

ro'' 



(14) 



(15) 



where vq is the lengthscale in the galactic plane. If we accept the usual form 
of softened isothermal distribution of neutrino dark matter |3^ 

where a is the core radius. Horizontal lengthscale is, then, given as 

'^o = (17) 

Thus, we obtain the general expression for flattening parameter in the "max- 
imal neutrino" halo assumption with additional ionizing sources: 



1 



na{anl — F)Tm^ 



(18) 



For a = 8 kpc, we obtain 



q = 6.601 < 



2.6 



Up 



(19) 



,10-2 cm-V 10-17 cm-3s-i 

Now we have cast the problem in the most tractable form. One may see that 
the following assumptions 

1. purely DDM ionization {F = 0), 

2. desirable lifetime of the unstable particle {T23 — 1), and 

3. Sciama's original value for the electron density of the Reynolds' layer 
(up = 0.04 cm-3), 
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truly lead to extreme flattening of the dark halo (g = 0.16). We would like to 
investigate the consequences of the relaxing of these assumptions. Keeping 
the assumption (1) and relaxing (2) and (3) leads to results shown in the 
Fig. 1. For the electron density, we use the more realistic values of the 
eq. i). 

We see that the "classical" value of T23 = 1 does not warrant an extreme 
degree of the neutrino halo flattening, contrary to the claim in the Ref. [p£ 
Instead, a moderate flattening q = 0.4 — 1 seems to be quite acceptable. 



5 Other ionizing sources 

If we relax the assumption (1), it is necessary to ask: which ionizing sources 
are viable for creating the electron density in (^? There are several possi- 
bilities, main ones being: 

• O and B stars of the galactic disk; 

• metagalactic UV ionizing background; 

• soft X-ray background, either of galactic or extragalactic origin; 

• cosmic rays which penetrate the Reynolds' layer; 

• large-scale ionizing shock waves. 

Some of them we can eliminate almost immediately, since we are not inter- 
ested in the detailed ionization structure, only in influence of a specific ioniz- 
ing source. Cosmic ray ionization, which has played such a prominent role in 
the first comprehensive ISM model of Field, Goldsmith and Habing in 1969 
|ri[] , is now considered to be of secondary importance: an average ionization 
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per H atom is now considered to be ^ 3 x 10^^^ s^^ (in contradistinction to 
earlier estimates of ~ 10~^^ s~^). This makes cosmic rays unimportant for 
the total ionization and heating budget of ISM, with important exception of 
the interiors of giant molecular clouds. If, as we have seen, the density of 
electrons is ~ 10^^ cm^^ (order of magnitude estimates), the total density 
can not be much higher, certainly < 10^^ cm^^. This makes resulting value 
of F in the eq. (0) quite unimportant in comparison with the DDM term, 
and resulting decrease in flattening is negligible. 

Similar considerations apply to soft X-ray ionization, which was first pro- 
posed by Silk and Werner The justification is that, like the cosmic rays. 
X-rays penetrate vast amount of matter before knocking out a photoelectron 
which is subsequently capable of ionizing and exciting hydrogen. Interest- 
ingly enough, for this ionization mechanism, more electrons are released from 
He than from H, since helium has larger photoionization cross-section at ener- 
gies ~ 100 eV, situation opposite from the one in case of DDM ionization. 
However, primary ionization rate for this process is estimated as about 10~^^ 
s~^ in an unshielded region Thus, resulting value of F is ^ 10^^^ cm~^ 
s~^, which does not significantly influence the required flattening of the neu- 
trino halo. The same conclusion applies a fortiori when we consider effects 
of finite opacity of the Reynolds' layer. 

What is the most difficult is to estimate the importance of the most 
plausible source of ionization: O and B stars in the Milky Way disk [|10| . This 
uncertainty is reflected in the conclusions of Mathis [|18|, and is connected 
with the (in) famous problem of the opacity of disks of spiral galaxies |1^ 

, which has not been solved to this day. It seems clear that the shape of the 
ionizing spectrum created by galactic early-type stars is capable of explaining 
the ionization state of the Reynolds' layer, but more detailed modelling of the 
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global ISM opacity will be necessary to positively discern whether this process 
really takes place. Our conclusion with respect to flattening seem, however, 
clear: if a fraction of O and B stars' ionizing flux successfully penetrates 
dense low-latitude ISM, resulting ionization would further obviate need for 
flattening of the hypothetical decaying neutrino halo. 

6 Plausible modifications of DDM theory 

Very recently, an attempt has been made to show that the diffuse ionization 
in the Galaxy can be explained via the decaying of sterile neutrinos with a 
mass of 27.4 eV and lifetime of ~ 10^^ s p^. Sterile neutrinos appear in the 



theory as a consequence of the recent Super-Kamiokande result according to 
which the atmospheric neutrino anomaly is mainly due to nearly maximal 
oscillations between and u^, where z/* is a sterile neutrino This ap- 



proach is conceived in order to overcome the perhaps biggest objection to the 
DDM theory, i.e. the bulk of the dark matter is assumed to be in neutrinos. 
Thus, according to Ref. |Q, one can state that matter contribution to the 



critical density of the Universe is Qm — 0.4, while the baryonic contribution 
is Qb — 0.08, that leads to halo density of sterile neutrinos of ~ 2 x 10^ 
cm~'^ that is 10 times less than the value previously assumed (cf. eq. 0]). 
According to the eq. ^ this would require that the radiative lifetime is sig- 
nificantly lowered for these sterile neutrinos, i.e. r ~ 10^^ s. However, as 
we have shown, it is not necessary to consider such radical change of this 
fundamental parameter of the DDM theory - it is possible to obtain lower 



concentrations of neutrinos by taking into account the term F, like in eq. [IC 
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7 Conclusions 



It is clear, from the eq. (|r^), that larger values for T23, require huge amount of 
flattening, which is quite unacceptable. For T23 ~ 10, corresponding flatten- 
ing parameter is g < 0.1. This much flattening will have several consequences 
contradictory to the observational data: it is extremely improbable that such 
number of microlensing events would have been observed (Ref. [0), it would 
violate the Oort limit (or its absence, see [0), causing either the excessive 
local halo density |T^, or too small total galactic mass to be consistent with 
the satellite galaxies measurements [^, ^ (but see [^). Both addi- 
tion of other ionizing sources and addition of other dynamically important 
halo components would only act to increase the discrepancy with the mini- 
mal tolerable amount of fiattening.|^ Thus, we consider this an indication of 
non-viability of the large T23 hypothesis. 

We conclude that in view of all the uncertainties considered, as well as 
improved values for several of the observable quantities, flattening of the dark 
halo of the Milky Way can be regarded as neither well-defined nor essential 
prediction of the DDM theory. 



The authors are happy to express their gratitude to Slobodan Ninkovic 
and Vesna Milosevic-Zdjelar for useful discussions. SS acknowledges the fi- 
nancial support of the University of Trieste. 



■^Strictly speaking, this will apply to the addition of another dynamically important 
halo component only if this additional component is not an ionizing source itself, which 
is certainly true for most of the envisaged components (CDM, MACHOs, cold gas, etc.). 
If we stretch our imagination and consider, for example, an additional species of decaying 
particles and call it X, then by choosing parameters such that T,x / (txttix) > ^i/fiT^m^) 
we could increase the resulting flattening parameter. However, any appeal of the original 
DDM (as well as conforming to the Occam's razor) would be lost in this case, which would 
pile unwarranted hypotheses and require fine-tuning of its model parameters. 
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